Određivanje uopštene kompleksne funkcije
assume(0 < R & 0 < C & 0 < L )
zamena = L == C*R^2
zamena = 
vrednosti = [C == 1, R == 1]
vrednosti = 
H(s) = simplify(subs(R/(R + 1/(1/(s*L) + s*C) + R), lhs(zamena), rhs(zamena)))
H(s) =

Amplitudska i fazna karakteristika
Hjw(w) = simplify(subs(H(s), s, 1i*w))
Hjw(w) =

Aw(w) = simplify(abs(Hjw(w)))
Aw(w) =

phi(w) = angle(Hjw(w))
phi(w) =

Crtanje amplitudske i fazne karakteristike
fplot(w, subs(Aw(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
fplot(w, subs(phi(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
Određivanje propusnog opsega
Aref = Aw(0)
Aref =

w3dB = expand(solve(Aw(w) == Aref/sqrt(2), w))
w3dB =

Aw_num(w) = subs(Aw, lhs(vrednosti), rhs(vrednosti))
Aw_num(w) =

wg1 = subs(w3dB(1), lhs(vrednosti), rhs(vrednosti))
wg1 =

wg2 = subs(w3dB(2), lhs(vrednosti), rhs(vrednosti))
wg2 =

Označavanje propusnog opsega u grafiku
fplot(w, subs(Aw(w), lhs(vrednosti), rhs(vrednosti)), [0,4])
rectangle('Position',[0, 0, double(wg1), double(Aref/sqrt(2))],...
'FaceColor',[1 0.9 0.9],...
rectangle('Position',[double(wg2), 0, double(4-wg2), double(Aref/sqrt(2))],...
'FaceColor',[1 0.9 0.9],...